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Abstract 

. We introduce notions of equivalence of conservation laws with respect to Lie symmetry 

groups for fixed systems of differential equations and with respect to equivalence groups or 
sets of admissible transformations for classes of such systems. We also revise the notion of 
linear dependence of conservation laws and define the notion of local dependence of potentials. 
To construct conservation laws, we develop and apply the most direct method which is 
effective to use in the case of two independent variables. Admitting possibility of dependence 
of conserved vectors on a number of potentials, we generalize the iteration procedure proposed 
by Bluman and Doran-Wu for finding nonlocal (potential) conservation laws. As an example, 
we completely classify potential conservation laws (including arbitrary order local ones) 
of diffusion-convection equations with respect to the equivalence group and construct an 
exhaustive list of locally inequivalent potential systems corresponding to these equations. 
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O ' 1 Introduction 

^t- ' 

After the Emmy Noether's remarkable paper [29] had become well-known, a number of authors 
■ (see e.g. [20,30,37,39]) searched for conservation laws using the symmetry approach based on 

the Noether's results. In view of the generalized Noether's theorem [30], there exists one-to-one 
i -^h | correspondence between the non-trivial generalized variational symmetries of some functional 

and the non-trivial conservation laws of the associated Euler-Lagrange equations, and any such 
symmetry is a generalized symmetry of the Euler-Lagrange equations. 

The Noether's approach has a number of advantages. It reduces construction of conservation 
laws to finding symmetries for which there exist a number of well-developed methods, and 
complete description of necessary symmetry properties is known for a lot of systems of differential 
equations. However, this approach can be applied only to Euler-Lagrange equations that form 
normal systems and admit symmetry groups satisfying an additional "variational" property of 
leaving the variational integral invariant in some sense [30]. The latter requirements lead to 
restriction of class of systems that could be investigated in such way. 

At the same time, the definition of conservation laws itself gives rise to a method of finding 
conservation laws. Technique of calculations used in the framework of this method is similar to 
the classical Lie method yielding symmetries of differential equations [19, Chapter 6]. As men- 
tioned in the above reference, such algorithmic possibility was first employed by P.-S. Laplace [26] 
for derivation of the well-known Laplace vector of the two-body Kepler problem. Following tradi- 
tion from group analysis of differential equations, we may call this method direct and distinguish 
four its versions, depending on the way of taking into account systems under investigation. (See 
e.g. [3,4,7,22] and Section 0] of this paper for more details as well as [43] for comparison of the 
versions and their realizations in computer algebra programs.) In the present paper we use the 
most direct version based on immediate solving of determining equations for conserved vectors 
of conservations laws on the solution manifolds of investigated systems. 
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Let us note that there exist other approaches for construction of conservation laws which 
differ from the Noether's or above direct ones, are based on exploitation of symmetry properties 
of differential equations and can be applied to non-Lagrangian systems. Thus, W.I. Fushchych 
and A.G. Nikitin [16] proposed to calculate directly bilinear combinations of solutions of motion 
equations, which are conserved in time by virtue of symmetries of these equations. It is possible 
in such way to find conservation laws corresponding to non-geometric symmetries. In the recent 
work [2] a purely algebraic formula has been derived for generating conservation laws of systems 
of differential equations that possess a scaling-invariance. 

To classify conservation laws, instead of the usual equivalence relation on their set (more 
exactly, on the set of conserved vectors) we propose to use the natural and more general notions 
of equivalences of conservation laws with respect to Lie symmetry groups for fixed systems 
of differential equations and with respect to equivalence groups or sets of admissible point 
(or contact) transformations for classes of such systems. Results of classification up to these 
equivalences are more comprehensible, especially, if a whole class of systems is studied and blend 
with the framework of group analysis. 

G. Bluman and P.R. Doran-Wu [7] proposed an ingenious procedure of branching iterations for 
finding nonlocal (potential) conservation laws of diffusion equations. Namely, on each iteration 
they use a conservation law from the previous iteration (one conservation law for one iteration) 
to introduce a potential and to construct the extended potential system. Then they study local 
conservation laws of the potential system, which are, generally speaking, nonlocal (potential) 
conservation laws for the initial equation. To the best of our knowledge, it was the first paper 
where the idea of hierarchy of potential systems and associated conservation laws is presented 
in an explicit form. 

We generalize the iteration procedure by G. Bluman and P. Doran-Wu, admitting dependence 
of conserved vectors on different number (from one to the maximum possible that) of new 
potentials on each iteration. The idea of a similar approach was adduced in [42] and was 
formalized in the form of notion of universal Abelian covering of differential equations [9,35,41]. 
Such approach naturally results in the questions on some independence of employed potentials. 
That is why, in this paper we consider definition of linear dependence of conservation laws in 
detail and define the notion of local dependence of potentials. 

Therefore, in the first part of the paper we propose some technique and discuss the classifica- 
tion problem for conservation laws in general. The ultimate goal of the second part is to present 
an exhaustive classification of potential conservation laws in a quite difficult and interesting 
case. As an illustration of the proposed technique, we choose the class of diffusion-convection 
equations 

ut = {A(u)u x ) x + B(u)u x , (1) 

where A = A(u) and B = B{u) are arbitrary smooth functions of u, A(u)^0. Symmetry prop- 
erties of were considered in [12,31,44], however the complete and strong group classification 
of Q was first presented in [33] (see also references therein for more details about symmetry 
analysis of classes intersecting class Q). 

Studying conservation laws of equations was started from linear equations [37]. 
V.A. Dorodnitsyn and S.R. Svirshchevskii [11] (see also [19, Chapter 10]) completely inves- 
tigated the local conservation laws for reaction-diffusion equations ut = (A(u)u x ) x + C(u). 
The first-order local conservation laws of equations Q were constructed by A.H. Kara and 
F.M. Mahomed [23]. Developing results of [7] obtained for the case B = 0, N.M. Ivanova [21] 
classified the first-order local conservation laws for equations (^Q) with respect to the equivalence 
group and constructed potential conservation laws of the first level. Namely, she made two steps 
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of the iteration procedure, looking, in the second step, for the first-order local conservation laws 
of the potential systems obtained after the first step. 

In the present work we exhaustively classify, with respect to the corresponding equivalence 
group, the local conservation laws of an arbitrary order, find the simplest and general poten- 
tial conservation laws and construct locally inequivalent potential systems of equations (JJ). 
All possible steps of the branching iteration procedure are done, and admission of dependence 
of conserved vectors on a number of potentials is of fundamental importance for completing 
iterations. We obtain eight inequivalent cases of equations Q having different sets of potential 
conservation laws: 

• the general case (the parameter- functions A and B are arbitrary); 

• three cases with an arbitrary value of A and a special value of B (B = 0; B = A] 
B = fAdu + uA); 

• three corresponding linearizable equations (A = u~ 2 , B = 0; A = B = u~ 2 ; A = 1, 
B = 2u) and 

• the linear heat equation (A = 1, B = 0). 

The latter case takes on special significance in our consideration since investigation of the 
linearizable equations having infinite series of potential conservation laws is reduced to this case 
and the non-linearizable equations from class (^Q) have at the most two independent conserva- 
tion laws. 

Our paper is organized as follows. First of all (Section EJ) we give a basic theoretical back- 
ground, following the spirit of [30]. We recall the notions of equivalence of conserved vectors 
and characteristics with respect to the triviality relation and discuss properties of the space of 
conserved vectors and the space of characteristics. This naturally results in the notions of linear 
dependence of conservation laws and local dependence of potentials (Sections 0] and EJ). 

In Section |3] we introduce the notions of equivalence of conservation laws with respect to 
Lie symmetry groups for fixed systems of differential equations and equivalence groups or sets 
of admissible point (or contact) transformations for classes of such systems. We emphasize 
possibility of solving classification problems for conservation laws with respect to the above 
equivalences similarly to usual group classification problems for differential equations. 

In Section 21 we adduce different versions of the direct method of construction of conservation 
laws, emphasizing the most direct one and combining it with classification up to symmetry or 
equivalence groups. We also generalize the iteration method of finding potential conservation 
laws and associated potential systems. 

Since the two-dimensional case is special for construction of potential systems we describe it 
in more detail in Sectional 

Then we apply the theoretical background given in the previous sections to investigation of 
diffusion-convection equations from class The local conservation laws of are classified 
with respect to the corresponding equivalence group in Section H3 In Section [7| we construct 
simplest potential conservation laws and analyze connections between them using potential 
equivalence transformations. The subject of Section |H1 is the description of general potential 
conservation laws of the linear heat equation. In Section |§] we complete studying potential 
conservation laws of diffusion-convection equations and adduce the hierarchy of conservation 
laws obtained in our framework. In the same section we give an exhaustive list of locally 
inequivalent potential systems of equations The obtained results can be interpreted as 

construction of universal Abelian covering for the whole class of diffusion-convection equations. 
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2 Basic definitions and statements 



Let £ be a system L(x,u^) = of I differential equations L 1 = 0, . . . , I) = for m unknown 
functions u = (u , . . . , ii m ) of n independent variables x = (x±, . . . , x n ). Here denotes the set 
of all the derivatives of the functions u with respect to x of order no greater than p, including 
u as the derivatives of the zero order. Let denote the set of all algebraically independent 
differential consequences that have, as differential equations, orders no greater than k. We 
identify with the manifold determined by in the jet space 

Definition 1. A conservation law of the system C is a divergence expression Div F := D^F 1 
which vanishes for all solutions of C: DivF|£= 0. The n-tuple F = (F 1 (x, U( r )), . . . , F n (x, U( r ))) 
is called a conserved vector of this (local) conservation law. 

In Definition ^ an d below Di = D Xi denotes the operator of total differentiation with respect 
to the variable Xj, i.e. Di = d x . + u a ai d u ^ where u a a and u a ai stand for the variables in jet 
spaces, which correspond to derivatives d^u a /dx^ 1 . . . dx" n and du%/dxi, a = («i,...,a n ), 
oii G N U {0}, |or|: = a\ + • ■ ■ + a n . We use the summation convention for repeated indices and 
assume any function as its zero-order derivative. The notation V^L means that values of V are 
considered only on solutions of the system C. 

Definition 2. A conserved vector F is called trivial if F l = F l + F l , i = l,n, where F l and 
F l are, likewise F l , functions of x and derivatives of u (i.e. differential functions), F l vanish 
on the solutions of C and the n-tuple F = (F , . . . ,F n ) is a null divergence (i.e. its divergence 
vanishes identically). 

The triviality concerning the vanishing conserved vectors on solutions of the system can be 
easily eliminated by confining on the manifold of the system, taking into account all its necessary 
differential consequences. 

A characterization of all null divergences is given by the following lemma (see e.g. [30]). 

Lemma 1. The n-tuple F = (F , . . . , F n ), n >2, is a null divergence (Div F = 0) iff there exist 
smooth functions v %:> = l,n) of x and derivatives of u, such that v 13 = —v 31 and F l = Djv 13 . 

The functions v %3 are called potentials corresponding to the null divergence F. If n = 1 any 
null divergence is constant. 

Definition 3. Two conserved vectors F and F' are called equivalent if the vector- function F' — F 
is a trivial conserved vector. 

Note 1. Sometimes other definitions of equivalence of conservation laws are used [15]. 

By the latter definition any trivial conserved vector is equivalent to the vanishing one. For any 
system C of differential equations the set CV(£) of conserved vectors of its conservation laws is 
a linear space, and the subset CVo(£) of trivial conserved vectors is a linear subspace in CV(£). 
The factor space CL(£) = CV(£)/ CVo(£) coincides with the set of equivalence classes of CV(£) 
with respect to the equivalence relation adduced in Definition |3J We can identify elements 
of CL(£) with conservation laws and call CL(£) also as the space of conservation laws of L (see 
e.g. [45]). That is why we assume description of the set of conservation laws as finding CL(£) 
that is equivalent to construction of either a basis if dimCL(£) < oo or a system of generatrices 
in the infinite dimensional case, and we will additionally identify elements from CL(£) with their 
representatives in CV(£). In contrast to the order rp of a conserved vector F as the maximal 
order of derivatives explicitly appearing in F, the order of a conservation law as an element T 
from CL(£) is called rnm{rp\F G CV(£) corresponds to J 7 }. Under linear dependence of 
conservation laws we understand linear dependence of corresponding elements in CL(£). 
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Definition 4. Conservation laws of a system £ are called linearly dependent if there exists their 
linear combination having a trivial conserved vector. 

Let the system £ be totally nondegenerate [30]. Then application of the Hadamard lemma 
to the definition of conservation law and integrating by parts imply that the left side of any 
conservation law of £ can be always presented up to the equivalence relation as a linear com- 
bination of left side of independent equations from £ with coefficients being functions on a 
suitable jet space 

jik). 

Div F = A M L M . (2) 
Here the order k is determined by £ and the allowable order of conservation laws, [i = 1,1. 

Definition 5. Formula © and the £-tuple A = (A 1 , . . . , A') are called the characteristic form 
and the characteristic of the conservation law Div F = correspondingly. 

The characteristic A is trivial if it vanishes for all solutions of £. Since £ is nondegenerate, 
the characteristics A and A satisfy (J2J) for the same F and, therefore, are called equivalent iff 
A — A is a trivial characteristic. Similarly to conserved vectors, the set Ch(£) of characteristics 
corresponding to conservation laws of the system £ is a linear space, and the subset Cho(£) of 
trivial characteristics is a linear subspace in Ch(£). The factor space Chf(£) = Ch(£)/ Cho(£) 
coincides with the set of equivalence classes of Ch(£) with respect to the above characteristic 
equivalence relation. 

The following result [30] forms the cornerstone for the methods of studying conservation laws, 
which are based on formula (j2j), including the Noether theorem and the direct method by Anco 
and Bluman [3,4]. 

Theorem 1 ( [30]). Let £ be a normal, totally nondegenerate system of differential equations. 
Then representation of conservation laws of £ in the characteristic form (J2J generates a one- 
to-one linear mapping between CL(£) and Chf(£). 

Using properties of total divergences, we can exclude the conserved vector F from (J2J) and 
obtain a condition for the characteristic A only. Namely, a differential function / is a total 
divergence, i.e. / = Div F for some n-tuple F of differential functions iff E(/) = 0. Hereafter 
the Euler operator E = (E 1 , . . . , E m ) is the m-tuple of differential operators 



E a = {—D) a d u a , a = l,m, 

where a = (ax, . . . ,a n ) runs the multi-indices set (a* <ENU{0}), (— D) a = (-Di) ai . . . (-D m ) am . 
Therefore, action of the Euler operator on (J2J) results to the equation 

E(A"L") = D* X (L) + Dl(X) = 0, (3) 

which is a necessary and sufficient condition on characteristics of conservation laws for the sys- 
tem £. The matrix differential operators D* x and are the adjoints of the Frechet derivatives D A 
and D^, i-e. 

= ((-»>■ (s^))- D " A » - (<-">" (S A " 

Since D^(£) = automatically on solutions of £ then equation (j3J) implies a necessary condition 
for A to belong to Ch(£): 

D2(A)| £ =0. (4) 

Condition Q can be considered as adjoint to the criteria C 7 ?) | ^ = fo r infinitesimal invariance 
of £ with respect to evolutionary vector field having the characteristic n = (rj , . . . ,rj m ). That 
is why solutions of ® are called sometimes as cosymmetries [6,35] or adjoint symmetries [4]. 
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3 Equivalence of conservation laws 



We can essentially simplify and order classification of conservation laws, taking into account 
additionally symmetry transformations of a system or equivalence transformations of a whole 
class of systems. Such problem is similar to one of group classification of differential equations. 

Proposition 1. Any point transformation g maps a class of equations in the conserved form 
into itself. More exactly, the transformation g: x = x g (x,u), u = u g (x,u) prolonged to the 
jet space transforms the equation DiF % = to the equation D^F* = 0. The transformed 
conserved vector F g is determined by the formula 

D x- 1 

F l g (x, = * J * F 3 (x,U( r )), i.e. F g (x,u^) = -—^-(D x x)F(x,u^) (5) 

in the matrix notions. Here |-t>a;x| is the determinant of the matrix D x x = (D Xj Xi). 

Note 2. In the case of one dependent variable (m = 1) g can be a contact transformation: 
x = x g (x,U(i}), = u g n\(x, Similar notes are also true for the statements below. 

Definition 6. Let G be a symmetry group of the system C. Two conservation laws with the 
conserved vectors F and F' are called G- equivalent if there exists a transformation g € G such 
that the conserved vectors F g and F' are equivalent in the sense of Definition |3J 

Any transformation g £ G induces a linear one-to-one mapping g* in CV(£), transforms 
trivial conserved vectors only to trivial ones (i.e. CVo(£) is invariant with respect to g*) and 
therefore induces a linear one-to-one mapping g^ in CL(£). It is obvious that gi preserves linear 
(in)dependence of elements in CL(£) and maps a basis (a set of generatrices) of CL(£) in a 
basis (a set of generatrices) of the same space. In such way we can consider the G-equivalence 
relation of conservation laws as well-determined on CL(£) and use it to classify conservation 
laws. 

Proposition 2. If the system C admits a one-parameter group of transformations then the 
infinitesimal generator X = £*<9j + rfd u °. of this group can be used for construction of new 
conservation laws from known ones. Namely, differentiating equation (J5J) with respect to the 
parameter e and taking the value e = 0, we obtain the new conserved vector 

F = -X (r) F* + (DjftF* - {D^)F\ (6) 

Here denotes the r-th prolongation [30,32] of the operator X. 

Note 3. Formula © can be directly extended to generalized symmetry operators (see, for 
example, [5,23]). A similar statement for generalized symmetry operators in evolutionary form 
(£ l = 0) was known earlier [18, 30]. It was used in [24] to introduce a notion of basis of 
conservation laws as a set which generates a whole set of conservation laws with action of 
generalized symmetry operators and operation of linear combination. 

Proposition 3. Any point transformation g between systems C and C induces a linear one- 
to-one mapping g* from CV(£) into CV(£), which maps CVo(£) into CVo(£) and generates a 
linear one-to-one mapping gi from CL(£) into CL(£). 

Corollary 1. Any "point transformation g between systems C and C induces a linear one-to-one 
mapping gf from Chf(£) into Chf(£). 
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It is possible to obtain an explicit formula for correspondence between characteristics of C 
and t. Let & = A^ v L v , where A» u = A^ va D a , A^ va are differential functions, a = (a x , . . . , a n ) 
runs the multi-indices set (aj €N U {0}), [/,, v = 1, I. Then 

Here A v ^* = (— D) a ■ A^ va is the adjoint to the operator A VfJi . For a number of cases, e.g. if C 
and C are single partial differential equations (1 = 1), the operators A^ u are simply differential 
functions (i.e. A^ va = for \a\ > 0) and, therefore, A"^* = A^. 

Consider the class C\s of systems L(x, ur p ), 9(x, utp))) = parameterized with the parameter- 
functions 9 = 9(x,U( p \). Here L is a tuple of fixed functions of x, U( p ) and 9. 9 denotes the 
tuple of arbitrary (parametric) functions 9(x,U( p \) = (9 1 (x,ui p \), . . . ,9 k (x,u^)) satisfying the 
condition S(x, U( p \,9r q -\(x, Uf p \)) = 0. This condition consists of differential equations on 9, where 
x and U(p) play the role of independent variables and Qr q \ stands for the set of all the partial 
derivatives of 9 of order no greater than q. In what follows we call the functions 9 arbitrary 
elements. Denote the local transformations group preserving the form of systems from C\s as 
G~ = G~(L,S). 

Consider the set P = P(L,S) of all pairs each of which consists of a system from C\s 
and a conservation law of this system. In view of Proposition 03 action of transformations 
from G~ together with the pure equivalence relation of conserved vectors naturally generates 
an equivalence relation on P. Classification of conservation laws with respect to G~ will be 
understood as classification in P with respect to the above equivalence relation. This problem 
can be investigated in the way that is similar to group classification in classes of systems of 
differential equations. Namely, we construct firstly the conservation laws that are defined for 
all values of the arbitrary elements. (The corresponding conserved vectors may depend on the 
arbitrary elements.) Then we classify, with respect to the equivalence group, arbitrary elements 
for each of that the system admits additional conservation laws. 

In an analogues way we also can introduce an equivalence relation on P generated by all 
admissible point or contact transformations (called also form-preserving ones [25]) in pairs of 
equations from C\s- 

Note 4. It can be easy shown that all the above equivalences are indeed equivalence relations. 

4 Direct iteration method of finding conservation laws 

To construct conservation laws of a system C of differential equations, we iterate a modification 
of the most direct method based on Definition More precisely, the algorithm is as follows. 

Zeroth iteration. At first we construct local conservation laws of C. We fix an (arbitrary) order r 
of conserved vectors under consideration. Then we introduce local coordinates ("unconstrained 
variables") on the manifold £( r +i) determined by the system C and its differential consequences 
in j( r+1 \ The other ("constrained") variables of j( r+1 ) are expressed via unconstrained ones by 
means of using the equations of Cr r +iy We substitute the expressions for constrained variables 
into a conservation law and split the obtained condition with respect to the unconstrained 
variables. This procedure results in a first-order linear system of determining equations for 
conserved vectors. Solving the determining equations up to the usual equivalence relation on 
CV(jC), we obtain complete description of local conservation laws of C. To classify conservation 
laws in easier and more systematic way (especially for classes of systems of differential equations), 
instead of usual equivalence we use the introduced above equivalence with respect to symmetry 
or equivalence transformations. (See Section El for examples.) 
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First iteration. After applying Lemma ^ to constructed conservation laws on the set of so- 
lutions of C = C°, we introduce potentials as additional dependent variables and attach the 
equations connecting the potentials with components of corresponding conserved vectors to C°. 
(If n > 2 the attached equations of such kind form an underdetermined system with respect to 
the potentials. Therefore, we can also attach gauge conditions on the potentials to C°.) 

We have to use linear independent conservation laws since otherwise the introduced potentials 
will be dependent in the following sense: there exists a linear combination of the potential tuples, 
which is, for some r' £ N, a tuple of functions of x and U( r n only. 

Then we exclude the unnecessary equations (i.e. the equations that are dependent on equa- 
tions of jC° and attached equations simultaneously) from the extended (potential) system C 1 
which will be called a potential system of the first level. Any conservation law of C° is a one 
of C 1 . We iterate the above procedure of the direct method for C 1 to find its conservation laws 
which are linear independent with ones from the previous iteration and will be called potential 
conservation laws of the first level. 

Further iterations. We make iterations while it is possible (i.e. the iteration procedure has to 
be stopped if all the conservation laws of a potential system C k+1 of the (k + l)-th level are 
linear dependent with the ones of C k ) or construct infinite chains of conservation laws by means 
of induction. This process may yield purely potential conservation laws of the initial system C, 
which are linear independent with local conservation laws and, therefore, depend explicitly on 
potential variables. 

Any conservation law from the previous step of iteration procedure will be a conservation 
law for the next step and vice versa, conservation laws which are obtained on the next step 
and depend only on variables of the previous step are linear dependent with conservation laws 
from the previous step. It is also obvious that the conservation laws used for construction of a 
potential system of the next level are trivial on the manifold of this system. 

Since gauge conditions on potentials can be chosen in many different ways, exhaustive real- 
ization of iterations is improbable in the case n > 2. 

The procedure of exclusion of constrained variables (which are described above in detail only 
for the zeroth iteration) is called in classical group analysis as "confining to the manifold of C . 
Taking into account C in the above way, we automatically eliminate the ambiguity connected 
with vanishing conserved vectors on the solutions of C However, the second kind ambiguity 
arising via existence of null divergences is preserved, and it is the main reason of difficulties in 
realization of this algorithm with symbolic computation systems [43] . 

To find conservation laws on each step of iteration procedure, one can apply other methods 
which are based on the characteristic form (J2J) or its consequences (J3J and (J3J) . These methods are 
also called as direct [3,4]. Following [43], for convenience we will numerate them as the second, 
third and fourth versions of the direct method in contrast to the above first one. They are 
close to the symmetry group method by Noether since in the case of Euler-Lagrange equations 
the coefficients A a are nothing else than Noether's characteristics. Taking into account the 
equivalence relation on Ch(£), one can assume during calculations that characteristics depend 
only on unconstrained variables. 

In the second version of the direct method the representation (J2J) is regarded as an equa- 
tion defined on an open subset of with respect to conserved vectors and characteristics 
simult aneously. 

In the framework of the third version, sought quantities are characteristics only. Determining 
equation Q is defined on an open subset of j( k \ Conserved vectors are reconstructed from 
known characteristics via explicit integral formulas. An algorithm of this (third) version of the 



8 



direct method was developed for Cauchy-Kovalevskaya systems by S. Anco and G. Bluman [3,4] 
(see also [30,45] for a theoretical background). 

The fourth version is based on equation @ which is defined on the manifold C and is only 
a necessary condition on characteristic of conservation laws. Therefore, one has to choose char- 
acteristics from the set of adjoint symmetries using additional conditions. Such approach was 
used by G. Bluman and P. Doran-Wu [7] for studying conservation laws of diffusion equations. 

Each from four above versions of the direct method has its advantages and disadvantages. 
A detailed comparative analysis of all the versions and their realizations in computer algebra 
programs are given by T. Wolf [43]. 

Prototypes and sources of a number of above ideas can be found in [27]. 

5 Two-dimensional case 

The case of two independent variables is singular, in particular, with respect to possible (con- 
stant) indeterminacy after introduction of potentials and high effectiveness of application of 
potential symmetries. That is why we consider some notions connected with conservation laws 
in this case separately. We denote independent variables as t (the time variable) and x (the 
space one). Any local conservation law has the form 

D t F(t,x,u (r) ) + D x G(t,x,u {r) ) = 0. (7) 

Here Dt and D x are the operators of the total differentiation with respect to t and x. F and G 
are called the conserved density and the flux of the conservation law correspondingly. 

Two conserved vectors (F,G) and (F',G') are equivalent if there exist such functions F, G 
and H of t, x and derivatives of u that F and G vanish on for some k and 

F' = F + F + D X H, G' = G + G- D t H. (8) 

Any conservation law of C allows us to deduce the new dependent (potential) variable v 
by means of the equations 

v x = F, v t = -G. (9) 

To construct a number of potentials in one step, we have to use a set of linear independent 
conservation laws (see the previous section) since otherwise the potentials will be dependent in 
the following sense: there exists a linear combination of the potentials, which is, for some r' € N, 
a function of t, x and Ui r i\ only. 

In the case of two independent variables we can also introduce the more general notion of 
potential dependence. 

Definition 7. The potentials v , ... ,v p are called locally dependent on the set of solution of the 
system C (or, briefly speaking, dependent) if there exist r' € N and a function H of the variables 
t, x, ur r i\, v , . . . , v p such that H{t,x,u^ r i^,v l , . . . , v p ) = for any solution (u, v , ... ,v p ) of the 
united system determining the set of potentials v , . . . , v p . 

Proof of local dependence or independence of potentials for general classes of differential 
equations is difficult since it is closely connected with precise description of possible structure 
of conservation laws. An example of such proof for diffusion-convection equations is presented 
below. 

In the case of single equation C, equations of form @ combine into the complete potential 
system since C is a differential consequence of 0- As a rule, systems of such kind admit a 
number of nontrivial symmetries and so they are of a great interest. 

Equations © and © imply the following statement. 



9 



Proposition 4. Any point transformation connecting two systems £ and £ of PDEs with two 
independent variables generates a one-to-one mapping between the sets of potential systems, 
which correspond to £ and £. Generation is made via trivial prolongation on the space of 
introduced potential variables, i.e. we can assume that the potentials are not transformed. 

Corollary 2. The Lie symmetry group of a system £ of differential equations generates an 
equivalence group on the set of potential systems corresponding to £. 

Corollary 3. Let C\$ be the set of all potential systems constructed for systems from the 
class £\s with their conservation laws. Action of transformations from G~(L,S) together with 
the equivalence relation of potentials naturally generates an equivalence relation on £\s- 

Note 5. Proposition^Jand its Corollaries imply that the equivalence group for a class of systems 
or the symmetry group for single system can be prolonged to potential variables for any step of 
the direct iteration procedure. It is natural the prolonged equivalence groups are used to classify 
possible conservation laws and potential systems in each iteration. Additional equivalences which 
exist in some subclasses of the class or arise after introducing potential variables can be used 
for deeper analysis of connections between conservation laws. 

6 Local conservation laws of diffusion— convection equations 

To classify the conservation laws of equations from class Q we have to start our investiga- 
tion from finding equivalence transformations. Application of the direct method to class ([T|) 
allows us to find the complete equivalence group G~ including the both continuous and discrete 
transformations. The following statement is true. 

Theorem 2 (see [33,34]). Any transformation from G~ has the form 

t = E^t + £\, X = £52 + Eft + E2, U = EqU + E3, A = E^ e\A, B = £4 1 £? ) B — £7, 

where £\, . . . , £7 are arbitrary constants, e^e^Eq ^ 0. 

The kernel of Lie symmetry group of equations from class (Q) is the group of the transforma- 
tions which are common for all of these equations. We denote it as G ker . Via trivial prolongation 
on the arbitrary elements A and B, G ker is isomorphic to a normal subgroup of G~. 

Theorem 3. G ker is formed by the transformations i = t + e\, x = x + £2, u = u, where £\ and 
£2 are arbitrary constants. 

First we search the conservation laws of equations from class Q in the form iff)). 

Lemma 2. Any local conservation law of any equation from class Q has the first order. More- 
over, up to equivalence on conserved vectors one can assume that the density depending on t, x 
and u and the flux depending on t, x, u and u x . 

Proof. Considering conservation laws on the manifold of equation Q and its differential conse- 
quences, we can assume that F and G depend only on t, x and Uk = d k u/dx k , k = 0,r', where 
r' < 2r. We expand the total derivatives in (J7J) and take into account differential consequences 
of the form u tj = D j x +2 JA + D x +1 J B, where / A = J A(u)du, fB = f B(u)du, j = 0~F. As a 
result we obtain the following condition 

F t + F u . (D x +2 J A + Di +1 JB) + G X + G U] u j+1 = 0. (10) 
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Let us decompose (fTU)) with respect to the highest derivatives Uj. So, the coefficients of u r / +2 
and iv+i gi ye the equations F u , = 0, G Ur , + AF U ,_ = that result in 

F = F, G = -Au rl F Urli + G, 

where F and G are functions of t, x, u, U\, . . . , u r /_i. Then, after selecting the terms containing 
u 2 r ,, we obtain that - AF Ur , _ ltt , _ x =0. It yields that F = F 1 u r >-x + F°, where F 1 and F° 
depend only on t, x, «, ui, . . . , u T >^2- 

Consider the conserved vector with the density F = F — D X H and the flux G = G + DfH, 
where H = f F 1 du r t_2- This conserved vector is equivalent to the initial one, and 

F = F(t, x, u, ui, . . . , u r '_2), G = G(t, x,u,ui, . . . ,u r >-i). 

Iterating this procedure a necessary number of times, we obtain the lemma statement. □ 

Note 6. A similar statement is true for an arbitrary (l+l)-dimensional evolution equation £ 
of the even order r = 2f, f £ N [1, 18]. For example [18], for any conservation law of C we can 
assume up to equivalence of conserved vectors that F and G depend only on t, x and derivatives 
of u with respect to x, and the maximal order of derivatives in F is not greater than f. 

Lemma El gives a stronger result for a more restricted class of equations. In the above proof 
we specially use the most direct method to demonstrate its effectiveness in quite general cases. 
This proof can be easily extended to other classes of (l+l)-dimensional evolution equations of 
odd orders and potential systems corresponding to equations from class (^Q) (see the proof of 
Lemma EJ). 

Theorem 4. Any equation from class Q has the conservation law Q where 



1. F = u, G = -Au x - JB. (11) 

A complete list of G~ -inequivalent equations Q having additional (i.e. linear independent 
with (|llj) ) conservation laws is exhausted by the following ones 

2. VA B = 0: F = xu, G = fA- xAu x , (12) 

3. VA B = A: F = (e x +e)u, G = —(e x + e)Au x — e J A, (13) 

4. A=l, B = 0: F = au, G = a x u-au x , (14) 



where e G {0,±1} mod G~, J A = J A{u)du, J B = jB(u)du, a = a(t,x) is an arbitrary 
solution of the linear heat equation at + a xx = 0. (Together with values A and B we also adduce 
complete lists of densities and the fluxes of additional conservation laws.) 

Proof. In view of Lemma|2 we can assume at once that F = F(t,x,u) and G = G(t,x,u,u x ). 
Let us substitute the expression for m deduced from Q into ® and decompose the obtained 
equation with respect to u xx . The coefficient of u xx gives the equation AF U + G Ux = 0, therefore 
G = —AF u u x + G l {t,x,u). Taking into account the latter expression for G and splitting the rest 
of equation (JJJ) with respect to the powers of u x , we obtain the system of PDEs on the functions 
F and G 1 of the form 

F uu = 0, BF U — AF UX + G\ = 0, F t + Gl = 0. (15) 
Solving first two equations of (fT5|) yields 

F = F 1 ^, x)u + F°(t, x), G 1 = F x fA - F 1 fB + G°(t, x). (16) 
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In further consideration the major role is played by the equation AF UXX — BF UX + F u t = that 
is a differential consequence of system (|15|) and can be rewritten as 




BFl+Fl = Q. 



Indeed, it is the unique classifying condition for this problem. There exist four different possi- 
bilities for values A and B. In all cases we obtain the equation F® + G x = 0. Therefore, up to 
conservation laws equivalence we can assume F° = G° = 0. Moreover, the function F 1 = const 
is a solution of the classifying condition in the general case. This solution corresponds to the 
conservation laws of Case 1. Only conservation laws of such kind exist for all admissible values of 
arbitrary elements A and B. Then we classify special values of A and B for which equation lfT|) 
possesses additional conservation laws. 

1. B (A, 1). Then = 0, F t = that gives contradiction with the assumption F 1 ^ const. 



2. A (1), B S (1). Then B = mod G~ and F xx = 0, F/ = 0, i.e. F 1 = x mod G~ 
(Case 2). 

3. B £ (A, 1),A,B# (1). Then B = A mod G~ and F xx + F x = 0, = 0, i.e. F 1 = e x + e 
mod G~, where e € {0,±1} (Case 3). 



Note 7. It follows from the proof that we can assume (A,B) ^ const in Cases 1, 2 and 3. (If 
(A, B) = const Cases 1, 2 and 3 are included in Case 4 for different values of a.) 

Using the conservation laws adduced in Theorem 0J we can introduce potentials for different 
values of the parameter-functions A and B and construct the corresponding potential systems 
(Cases 1-4 of Table 1). The important question for our consideration is whether the introduced 
potentials are locally independent in the sense of Definition If we know the precise structure 
of conservation laws the answer is almost obvious. 

Theorem 5. For any equation (^Q) potentials are locally dependent on the equation manifold iff 
the corresponding conservation laws are linear dependent. 

Proof. Since the direct statement of the theorem is obvious (see Section |5J), we prove only the 
inverse statement, using the rule of contraries. Suppose that potentials v°, . . . , v p introduced 
with (independent) conservation laws of Cases 1-4 are locally dependent. The vanishing p means 
local triviality of v as a potential, i.e. v° can be expressed in terms of local variables and the 
corresponding conservation law is trivial. That is why it is sufficient to investigate only the 
special cases when the number of independent conservation laws is greater than 1. Therefore, 
p = 1 if either B = or B = A and p € N/{0} for the linear heat equation. 

Without loss of generality we can assume that there exist r € N and a fixed function V 
of t, x, v = (v , . . . v p ) and U( r ) that v° = V(t,x,v,ur r \) for any solution of the united system 
determining the whole set of potentials v°, . . . , v p . Taking into account equation (J^) and its 
differential consequences, we can assume that V depends only on t, x, v and Uk = d k u/dx k , k = 
0, r', where r' < 2r. Let us apply the operator D x to the condition v° = V(t, x, v, u, ui, . . . , u r '): 
v x = V x + V v sv x + V u kUk+i- (Hereafter the index s takes the values from 1 to p.) Since in any case 
under consideration v x = f s u for some functions f s of t and x, we can split the differentiated 
condition with respect to Uk step-by-step in the reverse sequence, beginning with the major 
derivative. As a result, we obtain V u u = 0, V x = and f° = V v sf s , i.e. the functions /°, . . . f p 
are linear dependent. This gives a contradiction with the supposition that the conservation laws 
are independent. □ 



4. A,B e (1). Therefore, A = 1, B = mod G~ and F/ + F xx = (Case 4). 



□ 
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7 Simplest potential conservation laws 
of diffusion— convection equations 



Let us investigate local conservation laws of potential systems 1-4 from Table 1, which have the 
form 

D t F(t,x,u^,v^) + D x G(t,x,U( r ),V( r )) = 0. (17) 

These laws can be considered as nonlocal (potential) conservation laws of equations from class (^Q) . 
We assume them as simplest potential conservation laws since the corresponding potential sys- 
tems are constructed from the initial equation Q with one conservation law only 

We classify conservation laws up to the equivalence relation with respect to the transformation 
group Cr~ which is a result of the trivial prolongation of the group G~ from Theorem |2 to the 
space of the potential v. 

Lemma 3. Any conservation law of form (|17|) for each of systems 1-4 from Table 1 has the zero 
order, i.e. it is equivalent to a law with a conserved density F and a flux G that are independent 
on the (non-zero order) derivatives of u and v. 

Proof. Consider any from the systems 1-4. Taking it and its differential consequences into 
account, we can exclude dependence of F and G on the all (non-zero order) derivatives of v and 
the derivatives of u containing differentiation with respect to t. The remain part of the proof is 
completely similar to the one of Lemma EJ □ 

In an analogous to Theorem 0] and more cumbersome way we prove the following statement. 

Theorem 6. A complete set of G^ r -inequivalent conservation laws of form (|17jl for equations 
from class Q is presented in Table 1 with a double numeration of cases. 



Table 1. Conservation laws and potential systems of convection-diffusion equations. 



N 


A 


B 


F 


G 


Potential system 


1 


V 


V 


u 


-Au x -fB 


v x 


= u, v t = Au x + f B 


1.1 


V 





V 


-J A 


v x 


— u, w x — v, w t = / A 


1.2 


V 


A 


e x v 


-e x J A 


v x 


= u, w x = e x v, w t =e x JA 


1.3 


V 


[A + uA 


e v 


-e°fA 


v x 


= u, w x = e v , w t = e" JA 


1.4 


u~ 2 





a 


a v u~ x 


v x 


= U, W x = cr, Wt = — (T„U _1 


1.5 


u- 2 


u- 2 


ae x 


a v u~ 1 e x 


v x 


= it, w x = ae x , Wt = — a v u~ 1 e x 


1.6 


1 


2u 


V 

ae 


V V 

a x e — due 


Vx 


= it, w x = ae v , Wt — aue v — a x e" 


2 


V 





xu 


[A — xAu x 


v x 


= XU, Vt — xAu x — JA 


2.1 


V 





x~ 2 v 


-x- 1 [A 


v x 


= xu, w x — x~ 2 v, w t — x~ x J A 


3 


V 


A 


(e x + e)u 


~(e x +e)Au x -sJA 


v x 


= (e x + e)u, v t = (e x + e)Au x + eJA 


3.1 


V 


A 


e x 


e x 


V x 


e x 

= (e + e)u, w x = v, 
\e x + e) 2 


[e x +e) 2V 


e x + e ^ A 


Wt 


e x 

I A 

e x + e J 


4 


1 





OLU 


a x u — au x 


v x 


= OLU, Vt = OLU x — OL x U 


4.1 


1 





(£).• 




v x 
w x 


= OLU, 



Here e G {0, ±1}, \A = [A(u)du, fB = fB(u)du; a(t,x), (3{t,x) and (j(t,v) are arbitrary solutions of the linear 
heat equation (at + a xx = 0, fit + P xx = 0, at + <?w = 0). In Case 1.3 we assume f B = u f A for a conservation 
law to have the adduced form. 
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Note 8. To prove Theorem we use all independent differential consequences of correspondent 
potential systems. In Table 1 for the double numerated potential systems we omit equations 
containing vt since they are only differential consequences of equations of these systems. 

Let us analyze connections between conservation laws and ones between potential systems, 
which arise due to additional (including purely potential) equivalence transformations in some 
special cases. Below we assume A {l,u~ 2 } mod G~ as a general value of A. 

General case. Equation in the general case has the unique linear independent local con- 
servation law (Case 1) with the conserved vectors (jF , G ) = (u, — Au x ). All conservation laws 
of the corresponding potential system 



v 



Au x + jB, (18) 



are trivial, i.e. in our framework equation Q of the general form admits only trivial potential 
conservation laws. 

B = 0. Any equation of such form has two linear independent local conservation laws (Cases 1 
and 2) with the conserved vectors (F , G 1 ) = (u, —Au x ) and (F 2 , G 2 ) = (xu, J A — xAu x ), and 
any conservation law is G ker -equivalent to one of them. Using these conservation laws, we can 
introduce two potentials v 1 and v 2 , where 

v\ = u, v\ = Au x , (19) 
v 2 = xu, v 2 = xAu x - J A. (20) 

Equations l|19|) and (|20|) considered separately form two potential systems for equation JQ) with 
vanishing B in unknown functions (u, v 1 ) and (u, v 2 ) correspondingly. The third potential system 
is formed by equations lfT§|) and (|2*U)) simultaneously, and three functions u, v 1 and v 2 are assumed 
as unknown. Each from systems (|19j) and (120H has one linear independent local conservation law 
(Cases 1.1 and 2.1). These conservation laws with conserved vectors (F ,G ) = (v 1 , — J A) 
and (F 21 ,G 21 ) = (x~ 2 v 2 , — x _1 J A) are simplest potential conservation laws for equation (|T|) 
with vanishing B and allow us to introduce "second-level" potentials w 1 and w 2 . As a result, 
we obtain two potential systems of the next level: 

w l x = v\ w} = JA, (21) 

1 1 A. (22) 



2 2 —2 2 2 

Vt = XU, WZ = X V , WT = X 



At the same time, the simplest potential conservation laws are trivial on the solution manifold 
of the united system H19j) - (|2*U1) since 

F n = D x {xv l - v 2 ), G n = -D t {xv l - v 2 ), 
F 21 = D x {v l - x- V), G 21 = -D t {v l - x- V). 

Moreover w 1 = xv 1 — v 2 , w 2 = v 1 — x~ l v 2 , i.e. systems (|2*T|) . (|2*2*|) and (|T!l|) - lj2T)|) are locally 
equivalent. It implies that really system (fT9*]) - (|2T?]) is generated by only three independent 
equations. We can choose e.g. the equations 



v x = u, v x = xu, xvi ~ v t 



Although system ()21() formally belongs to the second level, it is the most convenient for further 
investigation since it has the simplest form. 
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B = A. This case is analyzed in the similar way to the previous one. Any equation with B = A 
has the two-dimensional space of local conservation laws. Up to G ker -equivalence, we have two 
possibilities for conserved vectors (Cases 1 and 3): 

(F 1 , G 1 ) = (u, -Au x - J A) and 

(F 3 , G 3 ) = ({e x + e)u, -(e x + e)Au x - e J A). 

Using these conservation laws, we can introduce two potentials v 1 and t> 3 , where 

v\ = u, v\ = Au x + J A, (23) 

vl = {f + £ )u, v$ = {e x + e)Au x + e J 'A. (24) 

Equations (|23|) and (|24() considered separately form two potential systems for equation JQ) with 
B = A in unknown functions (u, v 1 ) and (u, v 3 ) correspondingly. The third potential system is 
formed by equations (|23|) and (|24() simultaneously, and three functions u, v 1 and v 3 are assumed 
as unknown. Each from systems (|23|) and ()24() has one linear independent local conservation 
law (Cases 1.2 and 3.1). These conservation laws with conserved vectors 



(F 12 , G 12 ) = {e x v 1 , -e x J A) and 

(F 31 ,G 31 )= ( - &X —[A] 

y ' \(e x +e) 2 e x + e J J 



are simplest potential conservation laws for equation (^Q) with B = A and allow us to introduce 
"second- level" potentials w 1 and w 3 . As a result, we obtain two potential systems of the next 
level: 

v x = u, w x = e x v l , wj = e x jA, (25) 
v 3 x = (e x + e)u, w 3 x = {eX e +£)2 v 3 , w 3 = -^—jA. (26) 

At the same time, the simplest potential conservation laws are trivial on the solution manifold 
of the united system (|23|) -(|!? H) since 

F 12 = D x ((e x + ey - v 3 ), G 12 = -D t ((e x + e)^ 1 - v 3 ), 
?3l _ n I '„,! y3 ^ /^3l _ n ( ' l y3 



F =D x [v 1 - -— , G 3i = -D t v 

e x + £ ) \ e x + £ 



Moreover 



w 1 = (e x + £)v l - v 3 , w 3 = v 1 



3 



e x + £ ' 



i.e. systems (|25|). ((2*p]l and (|23| 1 -(|24 |1 are locally equivalent. It implies that really system 1)23(1 - 
PH) is generated by only three independent equations. We can choose e.g. the equations 



(e x + e)n, (e 1 + - v 3 = e x fA. 



Although system (|25|) formally belongs to the second level, it is the most convenient for further 
investigation since it has the simplest form. 
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B = J A + uA. The initial potential system in Case 1.3 is reduced to Case 1.2 by means of 
the hodograph transformation 

t = t, x = v, v = x, u = u~ l , A = u~ 2 A, (27) 

and the conservation law 1.3 is transformed to the local one of Case 3 where e = 0. The 
same transformation extended to the potential w as w = —w + ve x also reduces the potential 
system 1.3 to 1.2. 

Linearizable equations. It is well known [8,10,13,14,17,38,40] that equations (J"Q) G~-equi- 
valent to ones of Cases 1.4, 1.5, 1.6 are linearized by a nonlocal (so-called potential equiva- 
lence [28,34]) transformations to the linear heat equation. That is why these equations stand 
out against the other diffusion-convection equations with having an infinite number of linear 
independent purely potential conservation laws. 

The u~ 2 -diffusion equation u t = (u~ 2 u x ) x is reduced to the linear heat equation [8] by 
the hodograph transformation (|27|). More exactly, ()27|) is a local transformation between the 
corresponding potential systems v x = u, vt = u~ 2 u x and v x = u, vt = u x constructed by means 
of using the "common" conservation law (Case 1). The u _2 -diffusion equation has, as a subcase 
of the case B = 0, two linear independent local conservation laws with the conserved vectors 

F 1 = u, G 1 = -u- 2 u x , (28) 
F 2 = xu, G 2 = -xu~ 2 u x - vT 1 (29) 

(Cases 1 and 2 of Table 1) and the infinite series of potential conservation laws (Case 1.4) 
additionally. Under the action of hodograph transformation (|27|1 the conservation law with con- 
served vector (j28fl is transformed to the trivial one of the linear heat equation with the conserved 
vector (1,0). And vice versa, the conservation law of the linear heat equation corresponding to 
the value a = 1 (Case 4) is transformed by (|27|) to the trivial one of the u _2 -diffusion equation 
with the conserved vector (1,0). The conservation law with conserved vector 1)29)1 is trivial on 
the manifold of potential system constructed by means of 1)28)1 . is equivalent to the one from 
Case 1.4 with a = 1 and is transformed to Case 4.1, where a = 1 and (3 = x. Case 1.4 is reduced 
by ((2*7)) to Case 4, where a = a. 

Since the equation u± = (u~ 2 u x ) x + u~ 2 u x is reduced to the -u~ 2 -diffusion equation by the 
local transformation t = t, x = e x , u = e~ x u, its conservation laws are connected with ones of 
the linear heat equation in a way which is similar to the previous case. 

The potential systems v x = u, «f = u% + u 2 and v x = u, vt = u x constructed with the 
"common" conservation law for the Burgers equation ut = u xx + 2uu x and the linear heat 
equation ut = u xx are connected via the transformation 

t = t, x = x, u = ue" , v = e v . 

(Here the tilde variables correspond to the Burgers equation.) Let us note that really the 
famous Cole-Hopf transformation [10, 17] (first found in [14]) linearizes the Burgers equation 
to the "potential" heat equation vt = v xx [28,34]. The Burgers equation has the "common" 
local conservation law (Case 1) and the infinite series of simplest potential conservation laws 
(Case 1.6). The above transformation between the potential systems induces the one-to-one 
mapping preserving a between the infinite series 1.6 and the one 4 of the "potential" heat 
equation vt = v xx . Then, the conservation law of form 4 with the function a for the "potential" 
heat equation vt = v xx is equivalent to the one with the function a for the heat equation 
ut = u xx , where a = a x . The conservation law of Case 1 for the Burgers equation is trivial on 
the manifold of the corresponding potential system and is mapped to trivial one of the system 
v x = u, v t = u x . 
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Linear heat equation. The linear heat equation ut = u xx has an infinite dimensional space of 
local conservation laws [11], which is generated by conserved vectors of the form 

(F a , G a ) = (au, a x u — au x ), 

where a = a(t,x) is an arbitrary solution of the backward linear heat equation at + a xx = 0. 
Using a fixed conservation law of such form, we can introduce the potential v a , where 

= au, vf = au x — a x u. (30) 

System ()30|) has one infinite series of conservation laws (Case 4.1) with conserved vectors 

{F a P} G a P) =((t) V a ,-a(£) U-(£) V a ) , (31) 



a J x \ a J x \ a . 

where (3 = (3(t,x) is an arbitrary solution of the backward linear heat equation f3t + (5 XX = 0. 
These conservation laws are simplest potential ones for the linear heat equation and allow us to 
introduce "second- level" potentials w al3 . As a result, we obtain potential systems of the next 
level: 

wf=(?) v a , wf = a(?) u + (P) v«. (32) 



a /x \ a /x \ a /t 



Consider the system 

v x = au, = au x — a x u, v@ = (3u, ff = fiu x — j3 x u (33) 

that is the union of two potential systems of form (|3()[) corresponding to the local conservation 
laws with the conserved vectors (F a ,G a ) and (F^,G^). In a similar way to the previous cases 
we can state that the second-level potential conservation law with conserved vector (|3*T|> is trivial 
on the solution manifold of system Q33|) since 

F a/3 = D x (J-v a - vA and G aP = -D t (J-v a - . 

Moreover, systems (|3*2*|) and (|3^|) are connected via the local substitution 

a.0 @ a 3 
a 

It implies that really system ()33(l is generated by only three independent equations. We can 
choose e.g. the equations 

a B a @ a ( P\ 

v x = au, vz. = pu, —v t — Vf = a — u. 

a ' \ a I „ 



As a result of our analysis, we can formulate the following statement. 

Theorem 7. For any non-linearized equation (^Q) and the linear heat equation the potential 
systems of the second level, which are constructed by means of using single conservation law of 
the simplest potential systems, are equivalent to first-level potential systems obtained with pairs 
of conservation laws. 
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8 Potential conservation laws of linear heat equation 



With respect to G~-equivalence the linear heat equation is the unique linear equation in class (JJ) . 
Investigation of its general potential conservation laws plays the major role in classification 
of potential conservation laws for linearizable equations in class (^Q) and, therefore, for whole 
class (0). (The simplest potential conservation laws are studied in the previous section.) 

As proved in Theorem 0J the linear heat equation has the infinite series of local conservation 
laws. Fixing an arbitrary p £ N and choosing p linear independent solutions a = (a 1 , . . . , a p ) of 
the backward linear heat equation, we obtain p linear independent conservation laws with the 
conserved vectors (F s , G s ) = (a s u, a.%11 — a s u x ). (Hereafter s = l,p.) In view of Theorem [5] the 
potentials v = (y 1 , . . . , v p ) introduced with these conservation laws by the formulas 

v s x = a s u, = a s u x — a s x u (34) 

are independent in the sense of Definitional 

For the linear heat equation the complete set of first level potential conservation laws is 
indeed the union set of conservation laws of systems 1)34(1 corresponding to all possible values 
of p and p-tuples a. The following theorem is true. 

Theorem 8. Any local conservation law of system (|34[) is equivalent on the manifold of sys- 
tem (|34|) to a local conservation law of the linear heat equation. 

Corollary 4. For the linear heat equation potential conservation laws of any level are equivalent 
to local ones on the manifolds of the corresponding potential systems, and potentials of any level 
are locally expressed via local variables t, x, ur r \ (for some r ) and potentials of the first level only. 

We present the proof of Theorem El as the chain of simple and nice lemmas. 

Lemma 4. Any local conservation law of system (|34l) is equivalent to that with the conserved 
vector (Ku,K x u — Ku x ) where the function K = K(t,x,v) is determined by the system 

K t + K xx = 0, a s K xv s - a s x K v s = 0. (35) 

Proof. Consider a local conservation law of system (j34|) in the most general form, where the 
conserved vector is a vector- function of t, x and derivatives of the functions u and v s from 
the zero order up to some finite one. Taking into account system (|34[1 and its differential 
consequences, we can exclude dependence of the conserved vector on the all (non-zero order) 
derivatives of v s and the derivatives of u containing differentiation with respect to t. Similarly to 
Lemma |21 we can prove that the reduced conserved vector (F,G) does not depend on (non-zero 
order) derivatives of u and, moreover, F = F(t, x, v), G = -a s F v s (t, x, v)u + G°(t, x,v). The 
function F and G° satisfy the system 

a s a s 'F vSyS , = 0, a s G° v s = 2a s x F v s + a s F xv s, F t + G° x = 0. 

Let us pass on to the equivalent conserved vector (F, G), where F = F + D X H, G = G — DfH 
and H = H(t, x, v) is a solution of the equations H x = —F, H t = G. (The variables v s is 
assumed as parameters in the latter equations.) Then F = Ku, G = K x u — Ku x . The function 
K = a s H v s depends on t, x and v and satisfy system □ 

Lemma 5. Let the solutions a s = a s (t,x) and [3 s = (3 s (t,x) of the (backward) linear heat 
equation satisfy the additional condition a x f3 s — a s (5 x = 0. Then for any i,j G N 

af#-a!/a? = 0. (36) 
Hereafter the subscripts i and j denote the i-th andj-th order derivatives with respect to x. 
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Proof. We make the proof by means of mathematical induction with respect to the value i + j. 

Equation (|36|) is trivial for i + j = 0, coincides with the additional condition for i + j = 1 
and is obtained from this condition by means of differentiation with respect to x if i + j = 2. 

Let us suppose that the Lemma's statement is true if i + j = m — 1 and i + j = m and prove 
it for i + j = m + 1. Acting on equation (|36j) where i + jf = m — 1 with the operator c?t + <9. 
and taking into account the conditions af + a s xx = and /3| + (3 XX = 0, we obtain the equation 



XX 



a s i+1 (3j +1 - a s j+1 (3? +1 = 0. 

Therefore, the statement is true for i' + j' = m + 1, 1 < i',f < m (i' = i + 1, j' = j + 1). 
It remains to perform the proof in the case i' = m + 1, f = (or equivalently i' = 0, j' = m+ 1). 
For these values of and j' the statement is produced by subtracting the induced above equation 
a mPi ~ a i@m = from the results of differentiation of equation (|31)|) where i = m, j = with 
respect to x. □ 

Let W(ip l , . . . , denotes the Wronskian of the functions tp 1 , . . . , ip l with respect to the 
variable x, i.e. W^ 1 , . . . ,ip ) = det((p 3 i ) i j = \- 

Lemma 6. The solutions tp 1 = tp 1 (t, x) , . . . , tp l = ip l (t, x) of a linear evolution equation hp = 
are linear dependent iff W^ 1 , . . . , <p l ) = 0. 

Proof. Since the equation Lip = is linear and evolution the operator L is the sum of dt and 
linear differential operator with respect to x with the coefficients depending on t and x. If the 
functions cp 1 , . . . , ip are linear dependent then the equality W((p l , . . . , <p ) = is obvious. Let 
us prove the inverse statement. 

In the case I = 2 the condition W^ 1 ,^ 2 ) = implies p? = Cp 1 , where C is a smooth 
function of t. Acting on the latter equality with the operator L, we obtain Ctp 1 = 0, i.e. 
C = const or p 1 = 0. In any case the functions p 1 and p 2 are linear dependent. 

Suppose W(p 1 , . . . , tp ) = 0. Without loss of generality we can assume W(p 1 , . . . , p l ~ x ) ^ 0. 
(Otherwise we consider a less value of I.) Then p l = C k p k , where C k are smooth functions of 
t and the superscript k runs from 1 to I — 1. Action of the operator L on the latter equality 
results in the equation C k p k = that implies, in view of the condition W{p l , . . . , p 1 ^ 1 ) ^ 0, 
C k = const. It gives the Lemma's statement. □ 

Lemma 7. // a|/5| — a^/3f = for < i < j < p then W(a l , . . . , a p , /3 s ') = for any s' . 

Proof. Let M?'- denote the (p — l)-th order minor of W(a, f3 s '), which is obtained by means of 
deletion of s'-th and (p + l)-th columns corresponding to the functions a s ' and /3 s ' and i-th and 
j-th rows. Let us multiply the equation af/3 s — a s /3 s = by (— iy+i+ s +p+1 Mf- and convolve 
with respect to the indices i and j. In view of the Laplace theorem on determinant expansion 
we obtain 

Here the sign means that the function a s is substituted instead of the function a s ' and 
we have summation over the index s. The Lemma's statement easily follows from the latter 
equation since for any fixed s / s' we have W(a\ aS ^ aS i , /3 s ) = 0. □ 

Lemma 8. The general solution of system (|35|) can be presented in the form 

K = a s H vS + /3°, (37) 

where H is an arbitrary smooth function of v, (3° = (3 (t,x) is an arbitrary solution of the 
backward linear heat equation. 
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Proof. In view of Lemma |1] the functions a s and (3 s = K v s satisfy the conditions of Lemma |S] 
and, therefore, the ones of Lemma and the variables v are assumed as parameters. Since 
a s are linear independent it implies K v a = C as a s , where C as are smooth functions of the 
variables v only. Hereafter the indices s, a and ? run from 1 to p. The expressions for the 
cross derivatives K v a v ^ = C^a s = C^a s result in the equation C%i = which can be 
easily integrated: C as = P^a for some smooth function P s of the variables v. Substituting the 
expressions for C as in the equations on K and integrating, we obtain K = a s P s + /3°, where 
f3° = (3°(t,x) is a solution of the backward heat equation. The latter equality and the equation 
a s K xv s — a s x K v s = together imply the equation (a x a ( — a a a x )(P^cr — P£) = 0. Analogously 
to Lemma we can state for any i,j € N 

(of a} - aja?)(i* - P£) = 0. (38) 

Let M.?- c denote the (jp — 2)-th order minor of W(a), which is obtained by means of deletion 
of a'-th and ?'-th columns corresponding to the functions a a ' and a ? ' and i-ih and j-th rows. Let 
us multiply the equation 1)38(1 by (— l)*+-?+ cr + ? and convolve with respect to the indices i 

and j. In view of the Laplace theorem on determinant expansion we obtain 

W{a\ a ^ a< " = 0- (39) 

Here the sign "W means that the functions a a and a q is substituted instead of the function a a ' 
and a ? correspondingly and we have summation over the indices a and For any fixed a ^ a' 
and ? ^ we have W^(«| af a^. Q! <r' ^^q,?') = 0- Since Wia) ^ in view of linear independence 
of the functions a s , equation ()39|) implies — P a ' , = 0, i.e. P s = H v s for some smooth 
function H of v. □ 

In view of Lemma |S1 the conserved vector (Ku,K x u — Ku x ) from Lemma H has the form 
(f3°u + D X H, p° x u - P°u x - D t H). 

The proof of Theorem |H1 is completed. 



9 Potential conservation laws of 

nonlinear diffusion— convection equations 

More general potential conservation laws than simplest ones are admissible only if the investi- 
gated system has 

• either more than one linear independent local conservation laws (and, therefore, we can 
introduce a number of different potentials for the first iteration step) 

• or non-trivial simplest potential conservation laws. 

As shown in Section it is possible in class Q only for special values of the parameter func- 
tions A and B. In view of results of Section [7| for the linearizable equations and Theorem |HJ we 
can formulate the following statement. 

Theorem 9. For any linearizable equation from class Q all potential conservation laws of the 
second level are equivalent on the manifold of the corresponding potential systems to potential 
conservation laws of the first level. 

To investigate completely the potential conservation laws of equations from class ifTj). it 
remains to study the subclasses with B = and B = A, equations from which have two 
independent local conservation laws, and the subclass B = J A + uA reduced to B = A by 
means of potential equivalence transformations. 
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Theorem 10. All potential conservation laws of any equation from class Q with B = or 
B = A are trivial on the manifold of the united potential systems (|19j) . (|2U|) and (|23|) , (|24|) 

constructed with pairs of independent local conservation laws. 

Proof. Consider the united system (fTU|) . ((2*0]) (or ("""3*j) . (HU). (Below we write down the differ- 
ences of the second case with the first one in brackets.) Similarly to Lemma [31 we can assume 
without loss of generality that F = F(t, x, v , v 2 ) and G = G(t, x, u, v 1 , v 2 ). Let us split the equa- 
tion DtF + D X G = on the manifold determining by the united system. Integration of one from 
the obtained equations results in the following expression for the flux G: G = —(QF) fA + G°, 
where G° = G°(t, x, v , v 2 ) and Q = d v i + xd v 2 (Q = d v i + e x d v 2). The other equations form 
the system on the functions F and G°: 

Q 2 F = 0, QG° = 0, Ft + G° x = 0, (QF) X + F v 2 = ( (QF) X - F v i = ). 

Therefore, F = F 1 v 1 + F°, and F , F° and G° are functions of t, y = x and u> = xv 1 — v 2 
(lo = e x v 1 —v 2 ) for which F® = Fy, F^r = — G^, F® = —Gy. The latter system implies existence 
of such function H = H(t, y,uj) that 

F 1 = H w (F^e 1 '^), F° = H y , G° = -H t . 

Then, F = D X H, G = —DfH, i.e. the conservation law is trivial. □ 

As shown above, there exists the following chain of local transformations between potential 
systems: 1.3 < — ► 1.2 < — > (|23|) . (|24j) . i.e. system 1.3 is locally equivalent to system (|23|) . ("24"|) . 
In view of this fact and Theorem [TO] we obtain the following statement. 

Theorem 11. On the manifold of the potential system 1.3 all potential conservation laws of 
any equation from class Q with B = j A + uA are trivial. 

Summarizing the above results, we note that up to G~ -equivalence the hierarchy of conser- 
vation laws (including local ones) for diffusion-convection equations (["Q) has the form: 

• the "common" local conservation law (Case 1) for arbitrary values of the parameter func- 
tions A and B; 

• two independent local conservation laws if B = (Cases 1 and 2) or B = A (Cases 1 
and 3); 

• one "common" local conservation law (Case 1) and one simplest potential that (Case 1.3) 
if B = JA + uA; 

• the infinite series of local conservation laws (Case 4) for the linear heat equation; 

• one "common" local conservation law (Case 1) and the infinite series of simplest potential 
conservation laws (Case 1.6) for the Burgers equation; 

• two independent local conservation laws for the -u -2 -diffusion equation (Cases 1 and 2) and 
the equation ut = (u~ 2 u x ) x + u~ 2 u x (Cases 1 and 3) as subcases of B = and B = A and 
the infinite series of simplest potential conservation laws (Cases 1.4 and 1.5) additionally. 

Note 9. Above we did not consider in an explicit form action of transformations from Lie 
symmetry groups on conservation laws of corresponding equations or potential systems. For the 
majority of cases this action is quite trivial. For example, we use translations with respect to x 
to normalize the constant e in Case 3. In Case 2 the same translations result in adding the 
"common" conservation law to the special one of this case. 
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A non-obvious connection between independent conservation laws can be established only 
for A = u~ 4 / 3 , B = (Case 1 and Case 2) by means of the transformation t = t, x = 1 — x~ 1 , 
u = x 3 u from the Lie symmetry group of the corresponding equation. This fact was first 
discovered in [23] in the framework of the "operator" approach. It should be mentioned that 
the values A = u' 4/3 , B = give rise to the equation which distinguishes from non-linear 
diffusion-convection equations Q by singular Lie symmetry properties. 

The Lie symmetry group (?_ of the linear heat equation contains infinite dimensional nor- 
mal subgroup Q°_ formed by the linear superposition transformations u = u + f(t,x), where 
/ = f(t,x) is an arbitrary solution of the same equation. Up to the equivalence relation, trans- 
formations from Cp_ act identically on the set of conservations laws of Case 4. Action of the 
finite dimensional factor group Q-/Q_ on this set induces the analogous factor group Q+jQ\ on 
the set of solutions of the backward linear heat equation, which is varied over by the parameter- 
function a. 

The hierarchy of conservation laws generates the complete set of locally inequivalent potential 
systems for the class under consideration: 

• "common" potential system (|TH|) (Case 1); 

• additional simplest potential systems (|2Uj) (Case 2) or (|24j) (Case 3) if B = or B = A 
correspondingly; 

• second level potential systems (|2T|) (Case 1.1) and (|25|) (Case 1.2) (which are really equiv- 
alent to the united potential systems of the first level) if B = or B = A correspondingly; 

• system ()34|) with arbitrary number of locally independent potentials for the linear heat 
equation. 

Potential symmetries arising for equations from Cases 1 and 1.1 of Table 1 were studied 
by C. Sophocleous [36]. Complete investigation of the potential system (fTK|) (Case 1) with the 
symmetry point of view was carried out in [34]. 

10 Conclusion 

The notions and methods proposed in the paper are simple and effective tools for investigation 
of both local and pure potential conservation laws. They can be applied to a wide range of 
physically interesting systems of differential equations. At the same time, there exist a number 
of unresolved problems, in particular, on determining the number of necessary iterations for 
construction of an exhaustive list of independent potential conservation laws or on connections 
of our framework with Wahlquist-Estabrook prolongation structures [42]. 

The adduced results for diffusion-convection equations can be developed and generalized in 
a number of directions. So, studying different kinds of symmetries (Lie, nonclassical, general- 
ized ones) of constructed potential systems, we may obtain the corresponding kinds of potential 
symmetries (usual potential, nonclassical potential, generalized potential ones). Let us note that 
investigation of generalized symmetries is natural for potential systems, since potentials intro- 
duced with equivalent conservation laws are related, in general, via transformations depending 
on derivatives of local dependent variables. Analogously, local equivalence transformations be- 
tween potential systems constructed for different initial equations result in nonlocal (potential) 
equivalence transformations for the class under consideration. In such way it is possible to find 
new connections between well-studied diffusion-convection equations [34]. We believe that the 
same approach can be used for investigation of wider classes of differential equations, e.g. vari- 
able coefficient diffusion-convection equations. We also plan to study conservation laws of more 
general structure (e.g. ones with pseudopotentials). 
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